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We theoretically investigate the entangled states of an atomic ensemble that can be obtained via 
cavity-feedback, varying the atom-light coupling from weak to strong, and including a systematic 
treatment of decoherence. In the strong coupling regime for small atomic ensembles, the system 
is driven by cavity losses into a long-lived, highly-entangled many-body state that we characterize 
analytically. In the weak coupling regime for large ensembles, we find analytically the maximum 
spin squeezing that can be achieved by optimizing both the coupling and the atom number. This 
squeezing is fundamentally limited by spontaneous emission to a constant value, independent of the 
atom number. 
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Harnessing entanglement in many-body systems is of 
fundamental interest (TJ] and is the key requirement for 
quantum enhanced technologies, in particular quantum 
metrology ;2|. In this respect, many efforts have been 
devoted to prepare entangled states in atomic ensembles 
because of their high degree of coherence and their po¬ 
tential for precision measurement. Spin squeezed states 
as well as number states have been produced follow¬ 
ing methods based either on coherent evolution in the 
presence of a non-linearity in the atomic field @-0|, o r 
on quantum non-demolition measurement J(J-|8[. Among 
methods of the first kind, cavity feedback |5|, [9j is one of 
the most promising: it has already allowed for the cre¬ 
ation of highly squeezed states jl] and the effective non¬ 
linearity introduced by the atom-cavity coupling can be 
easily switched off, making it very attractive for metrol¬ 
ogy applications. 

In this Letter, we analyze the entangled states that 
can be produced by cavity feedback in different coupling 
regimes from weak to strong, and derive the ultimate 
limits of the metrology gain, extending the optimization 
of squeezing to unexplored domains of parameters values. 
After optimization of both the coupling strength and the 
atom number, we find a maximum squeezing limit that 
depends only on the atomic structure. 

Cavity feedback relies on the dispersive interaction be¬ 
tween one mode of an optical cavity and an ensemble of 
three level atoms, e.g. alkali atoms with a hyperfine split¬ 
ting in the ground state (see Figure 1). The atom-cavity 
system is characterized by the atom-cavity coupling g , 
the cavity linewidth (HWHM) ft, the atomic detuning 
A and the spontaneous emission rate T with A ^ T. 
The dynamics of entanglement is governed by the two di¬ 
mensionless quantities C = g 2 /(ftT) and 0o = 2g 2 /(ftA). 
The cooperativity C gives the ratio between the number 


of photons emitted in the cavity mode to spontaneously 
emitted photons as it can be seen by a Fermi golden rule 
argument , and a large C is favorable to entanglement 
because it minimizes the role of spontaneous emission. 
The parameter </>o represents the cavity frequency shift, 
normalized to the cavity linewidth, when a single atom 
changes its hyperfine state. In the regime </>o 1, pho¬ 

tons leaking from the cavity precisely measure the atom 
number difference between the two hyperfine states and 
therefore destroy coherence between them. One could 
expect that this may prevent the apparition of entangle¬ 
ment. However, this is not the case and we identify the 
condition to produce entanglement in this regime and 
characterize the produced states. They appear to have 
potential for metrology as signaled by their quantum 
Fisher information. In the regime <po 1, spin coher¬ 
ence can be maintained and our calculations confirm that 
this regime is optimal for producing spin-squeezed states. 
One important result is that the maximum squeezing is 
limited by the ratio of the excited state linewidth to the 
hyperfine splitting that should be as small as possible. 
As <f>o/C = 2T/A, the condition T/A 1 allows to 
maintain 0o small while maximizing C. 

We consider N atoms, with two (hyperfine) ground 
states |0) and |1) equally coupled with a constant g 
and opposite detunings ±A to an excited manifold \e) 
by a single cavity mode (see FigH^i). We introduce 
the collective spin operators S x + iS y = |l)(0|i, 

S z = \ |l)(l|i — |0)(0|i obtained by summing the 

effective spin 1/2 operators for each atom. The initial 
atomic state is a coherent spin state, each atom being 
in an even superposition of |0) and |1). A single off- 

resonant cavity photon shifts the energies of these levels 
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in opposite directions by an amount Through these 
opposite light shifts, the energy difference between lev- 
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FIG. 1. Principle of the cavity feedback scheme. Left: 
Three-level atoms are coupled to a cavity mode of frequency 
lJ c - We suppose that the optical transitions between the two 
ground states and the excited state equally couple to the cav¬ 
ity mode with a coupling constant g. Thus, tuning the cavity 
as shown in the figure results in equal and opposite light-shifts 
for the two ground states. The width of the excited state is 
T. Right: The evolution of the atomic state is induced by 
shining light onto the cavity at a frequency uj p which is de¬ 
tuned from (jj c by 5. The mean photon number in the cavity 
depends on S z the atom number difference between the two 
hyperfine states, via an atom-induced cavity detuning ncj)oS z 
with o. When 5 = k and (fioVN is small compared to the 
cavity linewidth, the atom-induced fluctuation of the photon 
number is proportional to <j>oS z and the dynamics can be well 
approximated by an effective \S Z model. 


els | 0 ) and | 1 ) depends on the cavity photon number c^c, 
that depends on its turn on the population difference S z 
as the atoms change the index of refraction in the cavity. 
Spin squeezing in this scheme occurs in the following way: 
the atomic quantum noise in S z , induces fluctuations of 
the cavity field intensity (as shown in Figjl]), which dur¬ 
ing the evolution are imprinted into the phases of each 
atom, thus correlating S y with the population imbalance 
S z [ 22 [. Assuming low saturation of the optical transi¬ 
tion g 2 (c f c)/A 2 C 1 , we eliminate the excited manifold 
\e) and describe each atom within the | 0 )-| 1 ) subspace. 
The unitary evolution is governed by 

Ho/h = (5 + K(j)oS z )c^c + - c ) ( 1 ) 


where c annihilates a photon of the cavity mode, rj is the 
cavity pumping rate, S = ( uj c — uj p ) is the empty cav¬ 
ity detuning, and we already introduced </>o = 2 g 2 / (ftA) 
that is also the single-photon atomic light shift properly 
normalized. Cavity losses and spontaneous emission in¬ 
cluding the possibility to scattering outside the | 0 ) — | 1 ) 
subspace are described by jump operators: d c — \[2 ftc 
and lia 
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di,xa = Vf\\X)(a \ iC ; T x = ^( 4 ) 

AA CLrrrj 


here cr, a' = 0,1 and I / 0,1 label the internal state, 
i labels the atom and a a ' a are amplitudes that depend 
on the atomic structure and field polarization (TBl]. The 
operator di^ a , di^xa refer to Rayleigh and Raman 
processes for the atom i. For any given eigenstate of the 
atomic operator S z with eigenvalue m G [—A/ 2 , A/ 2 ], 
the cavity field reaches in a time 1/ft a steady state that 
is a coherent state of amplitude a(m) 


fteff + i(S + Kfom) 

where ft e ff is the cavity linewidth in presence of the atoms 


fteff = ft 1 


4 (-^Ray 


■ r Ram + r x ) 


N «r^ 


(6) 


As we are interested in times t ft -1 , we shall neglect 
the transient effects and assume that the cavity field is 
in steady state from the beginning of the evolution j23[. 
In alkali atoms with a small hyperfine splitting in the 
excited state, choosing the states mp = 0 for | 0 ) and 
11), 7r-polarized light and a detuning close to half of the 
hyperfine energy splitting, one finds opposite light-shifts 
for the two states and no Raman processes coupling |0) 
and |1). In this first example we therefore restrict our 
analysis to Rayleigh processes that commute with S z . 
Raman spin-flipping processes and scattering to other 
states will be discussed later on in particular in relation 
to spin squeezing. Under these conditions, with an = a 22 
and CL(j(ji = 0 for cr 7 ^ cr', one can calculate the atomic 
density matrix p, expressed in the tensor product basis 
|e) = |ei, 62 ,..., ejsf) where e* = 0 , 1 refers to the internal 
state of the i-th atom 

<?lH?2) = T (a ( m?i )| a(mr2 ))l+2^+(iV-|l? 1 -e 2 ||)r Ray t 

x e mt[|a(m ?1 ) | 2 (<5//^+0 o m ri ) -\a(m ? 2 ) | 2 (<5//t+^ 0 m £ - 2 )] 
x e -(|a(m e -G| 2 + |a(m e -J| 2 )||e 1 -r 2 ||^^ ^ 

where ||e || = J2i=i n l e d anc i m e = ||F|| — A/2. The first 
line in Q represents decoherence due to loss of photons 
(through cavity losses and spontaneous emission) that 
are entangled with the atoms. The second line repre¬ 
sents the unitary evolution, whereas the third line is a 
second contribution of spontaneous emission that tends 
to kill all the off-diagonal elements of the density matrix 
by projecting single atoms into |0) or | 1 ). To explore the 
apparition of entanglement in the evolution starting from 
a coherent spin state |' 0 (O)) = [(| 0 ) + \l))/y/2\ N , we cal¬ 
culate the change in the system purity after tracing out 
one atom [171. Il8l| that we note PC (_purity change ): 

pc = Tin, 2 ,...Ap 2 ] - TV 2 ,...,jv[(Tnp) 2 ] (8) 


If all the atoms are correlated, tracing one of them can 
strongly influence the purity. Indeed, one can show that 
PC > 0 implies that the state is not separable. PC’s 






















3 


c 

1 10 100 1000 



FIG. 2. (Color online) Map of entangled states at different 
times and values of the phase shift 0o. The purity change is 
shown in color and the isolines of the squeezing parameter are 
solid black lines. Parameters : (5 = ft, N = 50, ( 77 /ft) 2 = 10 -2 , 
cm — CL 22 — 1.0, a aa t = 0 for a / cr', A/F = 500. 


maximum value is 1/2 obtained for a Schrodinger cat 
state. In a purely hamiltonian model H = xS'f, one has 
PC = |[1 — (cosxt) 2 ^ -1 )] 0- The main advantage of 
the quantity ((HJ) , is that it only requires the calculation 
of a trace and it can be computed even for relatively 
large atom numbers despite the large size of the Hilbert 
space. In Fig. [2] we show PC, as a function of time 
and of foVN that is the cavity detuning induced by the 
quantum fluctuations of S z (as A S z = VN/ 2). On the 
same plot we show the isolines for the spin squeezing 
parameter 0 


N A 2 S± 

\W 2 


(9) 


where A 2 S± is the minimal variance of the collective spin 
orthogonally to the mean spin direction and | (S) | is the 
mean spin length. Spin squeezed states £ 2 < 1 appear 
for coupling values such that 0o\0/V < 1, this conclusion 
illustrated for 50 atoms in Fig HI holds for much larger 
atom numbers (see Figj4]). 

Strong coupling regime - The purity change PC in Figj2] 
detects a larger and larger region of entangled states as 
0o VN increases, even after very short evolution times. 
To understand their nature, let us first consider the case 
without spontaneous emission. In this case decoherence 
is only due to cavity losses and one can use the Fock ba¬ 
sis | m) to express the density matrix whose off-diagonal 
elements decay as (a0m)|<a(ra / )) 1+2 'W If 0o 1, any 
state with m ^ 0 shifts the cavity out of resonance so 
that the cavity is practically dark, while the only distin¬ 
guished state is the twin-Fock state m = 0 which does 
not detune the cavity. Coherences between this state and 
all the others vanish in a time to — K \ a ^=o)\^ = 2ft/?7 2 
after which the initial coherent spin state |0) = ^Jpo \m = 
0 ) + x/(l Zr Po)\^ ± ) 0 is mapped onto the mixture: 


p = p 0 \m = 0 )(m = 0| + (1 -po)|f/)- L )(V’- L | (10) 


The subspace |V ; ~ L )('0~ L | of states m / 0 is preserved, 
both from decoherence and unitary evolution, for a time 
ti ^ w | a (m = i )|2 = (^/t 7 2 )0§/4 by which the cavity starts 
to distinguish the next Fock state. At longer times more 
and more coherences between different Fock states are 
killed until a complete mixture is reached. Note that 
ti/to = 0o/8 1 for strong coupling. Interestingly, 

the long-lived state m with partially removed coher¬ 
ences is highly entangled, its Fisher information scaling 
as Ip = /fN 3 / 2 for large N. In Figj3]we show the 

purity-change and the Fisher information as a function 
of time, for 0o = 14 and N = 10 atoms. Rayleigh spon¬ 
taneous emission and cavity losses are included. Fisher 
information and PC reach those of the state (usd around 
the time nt = 250 ~ 1.25ft£o and stay close to these val¬ 
ues until Kt = 2 x 10 4 ~ 2 nti indicating that our picture 
still holds in presence of spontaneous emission for small 
atomic ensembles. (^From the last row of Eq. © we 
see however that there are density matrix terms that are 
very sensitive to spontaneous emission and decay in a 
timescale 1/A/T, suggesting that the non-Gaussian state 
m is probably limited to small numbers of atoms. 



FIG. 3. (Color online) Purity-change (blue dashed line) and 
Fisher information (red solid line) as a function of time for 
= 14, N = 10 in presence of Rayleigh scattering. The 
horizontal lines give the analytical predictions for Ip and PC 
obtained from the state C0- Other parameters are as in 
Figia 

Spin-squeezed states - We now concentrate on spin 
squeezing and large atom numbers. In Fig0] we show 
the squeezing parameter optimized over time £min> as 
a function of 0 qVN for N = 10 5 , in a realistic config¬ 
uration for 87 Rb where we choose the clock transition 
|F = 1,771^ = 0) — |F = 2 ,mF — 0), 7r-polarized light 
on D2 line and a detuning close to half of the hyperfine 
energy splitting so that there are symmetric couplings 
an = &22 an d no Raman processes ai 2 = & 2 i = 0 |15j . 
The red solid curve, calculated from © includes cavity 
losses and Rayleigh spontaneous emission while the red 
dashed curve includes only cavity losses. We see that 
spontaneous emission is here important only for small 
values of (j)oVN- The blue dash-dotted curve is an effec¬ 
tive H = x& z model which we derive for (j)oVN <C 1 and 
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FIG. 4. (Color online) Spin squeezing optimized over time as 
a function of fa faN in a realistic configuration for 87 Rb (see 
text) with N = 10 5 , A/T = 563.39 and (r]/n) 2 = 10~ 2 and 
5 = k. Red solid line: full model © with cavity losses and 
Rayleigh jumps (an = a 22 = 0.702). Red dashed line: full 
model without spontaneous emission. Blue dash-dotted line: 
effective model in the regime fafaN <C 1 and K e s ~ k with 
an — < 3-22 — 0.702, ao,i — 0 * 1,0 — 0, and ax,i — flx,o — 0.497. 
Horizontal gray dotted line: £^ in in the large N limit dm) 


K eff ~ k, in which we can include all the loss processes 
and that we can solve analytically 0- For large enough 
coupling TV -1 / 10 < fafaN this models predicts a best 
squeezing limited by cavity losses as found in @0 

e mi n = l(3) 4/5 N~ 2/5 ; i min = ^|(3) 1/5 V- 3 / 5 . 

( 11 ) 

On the other hand, if fafaN > 1 the squeezing is lost in 
the full model because of non linearities that are not in 
the effective model that is here out of its limit of validity. 
Nonlinear effects come into play when the cavity detuning 
to due quantum fluctuations nfa faN exceeds the effective 
cavity linewidth ft e ff (the blue region in Fig. [l] exceeds 
size of the fringe). As we show now, the equality con¬ 
dition between theses two quantities allows to introduce 
a critical atom number N c that distinguishes between 
two different regimes for spin squeezing. In each of these 
regimes, N < N c and N > 7V C , an appropriate effective 
H = x$z m °del can be derived in some parameter range. 
Using the condition K(j)ofaN < K e s can be written as 


fafaN 



< 1 ; N c = 


4A/r 


ttRay + flRam + 


( 12 ) 


where we have introduced aR ay = 2a crjCr , aR am = 
\da'fa 2 / da,a and a x = Ex#o,i I a Xa| 2 /<W N c corre¬ 
sponds to such number of atoms that photon losses due 
to atom-scattering equal those due mirror transmission. 

For N <C N c which is the situation of Figj4j the condi¬ 
tion G2D gives fafaN < 1. For fafaN > 1 the detuning 
induced by the quantum noise becomes larger than the 
cavity linewidth giving rise to nonlinear effects destroying 
squeezing. 

For N N c we always have fafaN <C In this 

regime dominated by absorption, the cavity linewidth 


increases linearly with the atom number / k ~ 
Nfa/faN~ c , faster than the atoms induced cavity detun¬ 
ing oc ^oV^-By deriving a second effective H = x$z 
model, for n and N N c for cavity losses and 

Rayleigh jumps 0 , we find that (i) the best squeezing 
becomes independent of fa for large fafaN and, most 
importantly, (ii) the best squeezing has a non-zero limit 
for N —>> oo 


AT—s>oo f 2anT \ f q n 

Cmin e I —I (13) 

We show the onset of this new regime as N is increased 
in Fig|5]for A/T = 10. 



FIG. 5. Spin squeezing optimized over time as a function 
of fafaN , for A/T — 10 and ( 77 /ft ) 2 — 10 -2 . Solid lines: 
full model J7J) with cavity losses and Rayleigh scattering with 
an = a 22 = 0.702. Dot-dashed lines: analytical results in 
the regime ft e ff ~ ft and fafaN <C 1. Dotted lines: analytical 
results in the regime K e s ft and fafaN 1. Horizontal 
gray dotted line: £^ in in the large N limit m- From top to 
bottom: N = 10 3 , 10 4 , 10 5 , 10 6 , 10 7 . 


Conclusions We predicts that highly-entangled many- 
body states driven by cavity losses can be prepared by 
cavity feedback in the strong coupling regime for small 
samples, and a very large amount of spin-squeezing is 
reachable for large atom numbers in the weak coupling 
regime. The spin-squeezing limit we find (fl3l) is a very 
small value for alkali atoms suggesting that there is still 
room for improvement in the experimental achievements. 
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